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Abstract — In  this  paper,  low-complexity  joint  power  assign¬ 
ment  algorithms  are  developed  for  multi-source  multi-destination 
relay  networks  where  multiple  sources  share  a  common  relay 
that  forwards  all  received  signals  simultaneously  to  destinations. 
In  particular,  we  consider  the  following  power  optimization 
strategies:  (i)  Minimization  of  the  total  transmission  power 
of  the  sources  and  the  relay  under  the  constraint  that  the 
signal-to-interference-plus-noise  ratio  (SINK)  requirement  of  each 
source-destination  pair  is  satisfied,  and  (ii)  Maximization  of  the 
minimum  SINK  among  all  source-destination  pairs  subject  to  any 
given  total  power  budget.  Both  optimization  problems  involve  K 
power  variables,  where  K  is  the  number  of  source-destination 
pairs  in  the  network,  and  an  exhaustive  search  is  prohibitive  for 
large  K .  In  this  work,  we  develop  a  methodology  that  allows  us 
to  obtain  an  asymptotically  tight  approximation  of  the  SINK  and 
reformulate  the  original  optimization  problems  to  single-variable 
optimization  problems,  which  can  be  easily  solved  by  numerical 
search  of  the  single  variable.  Then,  the  corresponding  optimal 
transmission  power  at  each  source  and  relay  can  be  calculated 
directly.  The  proposed  optimization  schemes  are  scalable  and 
lead  to  power  assignment  algorithms  that  exhibit  the  same  opti¬ 
mization  complexity  for  any  number  (  K)  of  source-destination 
pairs  in  the  network.  Moreover,  we  apply  the  methodology  that 
we  developed  to  solve  a  related  max-min  SINK  based  optimiza¬ 
tion  problem  in  which  we  determine  power  assignment  for  the 
sources  and  the  relay  to  maximize  the  minimum  SINK  among  all 
source-destination  pairs  subject  to  any  given  total  power  budget. 
Extensive  numerical  studies  illustrate  and  validate  our  theoretical 
developments. 

Index  Terms — Cooperative  networks,  interference  relay  chan¬ 
nels,  max-min  SINK  optimization,  multi-source  multi-destination 
relay  networks,  optimum  power  allocation. 


I.  Introduction 

RECENT  work  on  information-theoretic  aspects  of  coop¬ 
erative  relaying,  as  well  as  recent  proposals  of  practical 
cooperative  relaying  protocols  ([1]— [10]  and  references  therein) 
suggest  that  cooperative  relaying  may  lead  to  significant  im- 
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provements  in  detection  reliability  at  destinations  and  overall 
system  performance.  In  cooperative  relaying,  a  user/node  may 
serve  as  a  relay  and  assist  others  by  forwarding  their  signals  to 
destinations,  thus  enhancing  detection  reliability  at  the  destina¬ 
tion.  Various  relaying  strategies  have  been  studied  in  the  liter¬ 
ature  for  relay  to  forward  signals.  For  example,  relay  may  de¬ 
code  the  received  signal  and  forward  the  decoded  information  to 
a  destination,  or  it  may  simply  amplify  the  received  signal  and 
forward  it  to  the  destination. 

More  recently,  there  is  increasing  interest  in  investigating 
the  advantages  of  relaying  in  multi- source  multi-destination 
networks  [1 1] — [18],  [24]— [26],  which  promise  significant 
achievable  rate  improvement  in  shared- spectrum  multiple 
access  wireless  networks.  The  simplest  multi-source  multi-des¬ 
tination  relay  network  is  modeled  as  an  interference  relay 
channel  (IRC)  [11]  where  a  relay  helps  two  independent 
source-destination  pairs  by  using  different  relaying  strate¬ 
gies  such  as  decode-and-forward,  amplify-and-forward, 
or  compress-and-forward.  Past  literature  on  multi-source 
multi-destination  relay  networks  focused  primarily  on  infor¬ 
mation-theoretic  studies  including  achievable  rate  regions  or 
bounds  of  capacity  region  [1 1]— [18].  For  example,  in  [1 1]  a  rate 
splitting  technique  is  used  to  study  the  problem  of  achievable 
rate  region  for  a  Gaussian  IRC  channel,  where  each  message 
is  split  into  a  common  message  which  is  decodable  at  all  des¬ 
tinations  and  a  private  message  which  is  decodable  only  at  the 
intended  destination  [14].  Since  the  receivers  are  able  to  decode 
part  of  the  interference  messages,  the  effect  of  interference 
is  reduced  and  the  overall  communication  rate  is  therefore 
increased.  The  achievable  rate  region  of  [11]  was  further 
improved  in  [12]  by  considering  both  intended  message  and 
interference  forwarding  at  the  relay  (optimal  relay  strategies 
were  studied  under  the  assumption  that  the  relay  is  connected 
to  each  source  and  each  destination  via  orthogonal  and  finite 
capacity  links).  The  capacity  region  of  the  interference  channel 
with  a  single-relay  was  investigated  in  [15]  and  [16],  where  it 
was  shown  that  forwarding  the  intended  message  of  just  one 
source,  the  achievable  rates  for  both  source-destination  pairs 
can  be  improved.  By  assuming  that  the  relay  knows  the  source 
message  a  priori ,  a  relaying  strategy  was  proposed  in  [17]  and 
[18]  where  generalized  beamforming  with  dirty  paper  coding 
was  considered  for  a  two- source  two-destination  relay  network 
to  further  improve  the  capacity  region  of  the  network. 

Power  control  is  important  to  improve  overall  performance  in 
multi-source  multi-destination  relay  networks  [19].  While  past 
literature  offers  a  significant  amount  of  work  on  power  alloca¬ 
tion  for  single-source  single-destination  relay  networks  (see,  for 
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example  [20] -[23]  and  references  therein),  there  are  rather  lim¬ 
ited  studies  on  power  optimization  for  multi-source  multi-desti¬ 
nation  relay  networks.  In  [24],  power  allocation  was  optimized 
by  exhaustive  search  for  a  two- source  two-destination  relay  net¬ 
work  where  a  half-duplex  decode- and-forward  relay  was  con¬ 
sidered.  Unfortunately  the  exhaustive  search  is  not  scalable  and 
leads  to  prohibitive  optimization  complexity  for  networks  with 
larger  number  of  source-destination  pairs.  In  [25]  and  [26],  a 
power  allocation  scheme  was  proposed  for  a  multi-source  multi¬ 
destination  relay  network  based  on  geometric  programming  (the 
scheme  assumes  that  signals  from  different  sources  are  sent 
through  orthogonal  channels  and  the  direct  transmission  link  is 
not  involved  in  detection  at  destinations). 

In  this  paper,  we  analyze  and  optimize  a  general  multi-source 
multi-destination  relay  network  with  K  source-destination 
pairs  ( K  can  be  large).  The  network  allows  simultaneous 
multi-source  transmissions  through  nonorthogonal,  in  general, 
channels.  We  examine  two  power  optimization  strategies:  i) 
Minimization  of  the  total  power  consumption  of  all  sources 
and  relay  under  the  constraint  that  the  signal-to-interfer- 
ence-plus-noise  ratio  (SINR)  requirement  of  each  source-des¬ 
tination  pair  is  satisfied;  and  ii)  maximization  of  the  minimum 
SINR  among  all  source-destination  pairs  subject  to  any  given 
total  power  budget.  Thanks  to  an  asymptotically  tight  approxi¬ 
mation  of  the  SINR  that  we  develop,  we  are  able  to  reformulate 
the  original  optimization  problems,  which  involve  K  power 
variables,  to  single-variable  optimization  problems.  Then,  the 
resulting  optimization  problems  can  be  easily  solved  by  a 
simple  numerical  search  of  the  single  variable.  The  proposed 
optimization  schemes  are  scalable  and  lead  to  power  assign¬ 
ment  algorithms  that  exhibit  the  same  optimization  complexity 
for  any  number  (K)  of  source-destination  pairs  in  the  network. 
Moreover,  for  the  special  case  of  transmission  over  orthogonal 
channels,  we  are  able  to  further  simplify  the  single- variable 
optimizations  and  obtain  analytical  solutions  for  a  symmetric 
system.  Extensive  numerical  studies  included  in  this  paper 
illustrate  and  validate  our  theoretical  developments,  and  show 
that,  the  proposed  power  assignment  is  almost  identical  to  the 
exhaustive  search  method,  and  the  optimum  power  assignment 
schemes,  in  general,  can  significantly  improve  the  performance 
of  multi-source  multi-destination  relay  networks  compared  to 
an  equal  power  assignment  scheme. 

The  paper  is  organized  as  follows.  In  Section  II,  we  introduce 
briefly  the  system  model  of  a  multi-source  multi-destination 
relay  network  where  transmissions  occur  over  nonorthogonal, 
in  general,  channels.  In  Section  III,  we  determine  the  maximum 
ratio  combining  of  the  received  signals  at  each  intended  destina¬ 
tion  and  exploit  the  resulting  SINR.  In  Section  IV,  we  determine 
the  optimum  power  assignment  for  the  sources  and  the  relay  that 
minimizes  the  total  power  consumption  under  the  condition  that 
the  SINR  requirement  of  each  source-destination  pair  is  satis¬ 
fied.  In  Section  V,  we  determine  the  optimum  power  assignment 
that  maximizes  the  minimum  SINR  among  all  source-destina¬ 
tion  pairs  subject  to  any  given  total  power  budget.  Numerical 
studies  are  provided  in  Section  VI,  and  finally,  some  conclu¬ 
sions  are  drawn  in  Section  VII. 

The  following  notation  is  used  in  the  paper.  Bold  letters  in 
uppercase  and  lowercase  denote  matrices  and  vectors,  respec- 


Fig.  1.  Multi-source  multi-destination  relay  network. 


tively.  ( •  )*,  ( •  )T  and  ( •  )H  represent  the  conjugate,  the  trans¬ 
pose  and  the  Hermitian  transpose  operation,  respectively.  |  •  | 
and  ||  •  ||  represent  the  absolute  value  of  a  complex  number  and 
the  Frobenius  norm  of  a  vector/matrix,  respectively.  1^  is  an 
L  x  L  identity  matrix.  diag(/ii,  /12,  •  •  • ,  hi)  is  an  L  x  L  diag¬ 
onal  matrix  with  diagonal  elements  hi,  h^,  . . . ,  h^.  denotes 
a  sub-matrix  of  A  obtained  by  deleting  the  kth  column  and  kth 
row  of  A.  If  8ik  represents  the  kth  column  of  the  matrix  A,  then 
a*.  denotes  the  vector  obtained  after  removing  the  kth  entry  from 
3-k  • 

II.  System  Model 

For  illustration  purposes  and  simplicity  in  presentation,  we 
consider  a  single  relay  code  division  multiplexing  system  with 
K  sources  and  K  destinations  as  shown  in  Fig.  1,  where  trans¬ 
missions  occur  over  nonorthogonal,  in  general,  channels.  Our 
developments  can  be  generalized  to  multiple-relay  systems  and 
other  multiplexing  schemes  in  frequency  and/or  time.  Fet  Sk 
denote  the  kth  source  and  D\~  the  corresponding  destination, 
k  —  1, 2, . . . ,  K,  and  let  R  denote  the  relay.  The  relay  forwards 
simultaneously  the  signals  received  from  all  sources.  Fet  bk  de¬ 
note  the  transmitted  information  symbol  of  the  source  Sk  with 
unity  average  energy,  i.e.,  E{\bk\2}  =  1,  Vfc.  The  signal  sent  by 
the  source  Sk  can  be  expressed  as 

Sfc  —  C  kbki 

where  c&  =  (c^\  c®  •  •  • ,  c^)T  is  the  code/signature  of  the 
source  Sk ,  which  is  a  unit-energy  column  vector  with  length 
L.  The  codes/channels  of  different  sources  are,  in  general,  cor¬ 
related.  Fet  pkj  =  cl c j  denote  the  cross-correlation  between 
codes/channels  k  and  j,  where  pkj  £  [0, 1)  for  k  /  j,  and 

pkk  =  1.  Let  R  =  (pkj)  denote  the  K  x  K  cross-correlation 
matrix,  i.e. 

R=  (ci,c2,...,ck)T(ci,c2,...,ck).  (1) 

We  consider  the  following  two-phase  amplify-and-forward 
relay  strategy  with  L  time  slots  in  each  phase.  In  Phase  1,  each 
source  Sk  transmits  the  signal  with  transmitted  power  Pk. 
Then,  the  received  signals  at  the  destination  Dk  and  at  the  relay 
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R  during  the  7th  (1  <  i  <  L)  time  slot  can  be  modeled,  respec¬ 
tively,  as 

K 


y. 


1  =  1 
K 


/(*) 


=  E  VPlh^Jt 


(i)bi  +  nW,  l<i<L. 


(3) 


2=1 


In  Phase  2,  the  relay  amplifies  the  received  signals  and  forwards 
them  to  the  destination  with  an  amplification  factor  a  and  trans¬ 
mission  power  Pr.  The  received  signal  at  the  destination  Dk 
during  the  7th  (1  <  i  <  L)  time  slot  can  be  written  as 


VrL  =  VPrh^day^+n(X,  i=l,2,...,L.  (4) 

In  (2)-(4),  h^dk ,  h^r,  and  h^dk ,  (l,k  =  1, . . . ,  K )  are  the  co- 
efficients  of  the  channels  between  the  source  Si  and  the  desti¬ 
nation  Dk ,  between  the  source  Si  and  the  relay  R ,  and  between 
the  relay  R  and  the  destination  Dk ,  respectively,  during  the  7th 
(1  <  7  <  L)  time  slot.  n^\  and  represent  noise  at  the 
destination  Dk  during  the  7th  time  slot  of  Phase  1  and  Phase 
2,  correspondingly,  while  n^l  represents  noise  at  the  relay  R 
during  the  7th  time  slot.  The  channels  d  ,h^r,  and 
are  assumed  to  be  independent  Gaussian  random  variables  with 
zero-mean  and  variances  a2Sj  dk ,  (r2ir,  and  (T2  dk  „  respectively. 
All  noise  terms  are  assumed  to  be  independent  Gaussian  random 
variables  with  zero-mean  and  variances  a2 .  Without  loss  of  gen¬ 
erality,  we  assume  a2  =  1. 

The  channel  coefficients  in  matrix  format  can  be 
written  as  =  diag^J^,...,^^),  H„,r  = 

diag (h^r,...,hif}r)  and  Ur4k  =  diag(/i^,...,/i^J. 
Then,  the  received  signals  can  be  expressed  as  follows: 


k  _ 

Ys,dk  ~  Ev^Hai  ,dkcibi  +  n  (5) 

i=i 

K 

Ys,r  =  ^  ^  ,r^lbl  +  (6) 

1=1 

Yr,dk  =  \j Pr®P^-r,dkY  s,r  ^r,dk  i  (7) 


where  ys4k  =  y^t . . . ,  </^)T 

(ys]r ,  2/$ ,  •  •  • ,  yiLr  )T  and  yr  A  =  (yfy ,  y ^ , 

In  (5)-(7),  the  noise  vectors  ns?c/fc,ns?r  and  nr?^fc 
elements  that  are  independent  Gaussian  random  variables  with 
zero  mean  and  unit  variance.  The  amplification  factor  in  (7) 
is  specified  as 


y  s,r  = 

77(L)  )T 

yr,dk)  * 

have 


a 


2 


1 

S/= 1  PlPsi,r  +  ^ 


(8) 


where  /?Sz?r  =  cjT  E{ H  J[;T.HS/  ,r  }c/ .  By  substituting  (6)  and  (8) 
into  (7),  we  obtain 

K 

Yr,dk  =  V^^Hrv4  ^  ^/Pl'H.s^rClbl 
i=i 

H-  y/Pr<xHr^  n.s,r  +  nr?dfc .  (9) 


jointly  detect  the  information  transmitted  by  the  source  Sk-  The 
combined  signal  from  Phase  1  and  Phase  2  at  destination  Dk 
can  be  expressed  in  vector  form  as  follows: 


Yk 

where 


A  (  Ys,dk  \  _ 
Yr,dk 


H/  fc  = 


K 


H  k,k^kbk  +  H  i^ibi  +  n&, 

i=\d^k 


V^H  s^dk 
y/ P[PrO t  H  Sl  ?rH  r,dk 


is  a  2 L  x  L  virtual  channel  matrix  from  the  source  Sj  to  the 
destination  Dk ,  and 


nk  = 


y/  P raU7,dk  ns?r  -f  wr^dk 


is  an  equivalent  noise  vector  of  length  2 L.  We  note  that  H*.*.  is 
the  channel  matrix  associated  with  the  desired  source  Sk,  while 
H ik,  (/  /  k)  are  the  channel  matrices  of  the  interfering  sources. 
Based  on  maximum  ratio  combining  (MRC)  detection  [28],  the 
transmitted  signal  from  the  source  Sk  is  detected  as 

bk  =  arg  min  |w?0yfc  —  bk\2 ,  (10) 

bk  £A 


where  A  is  the  set  of  transmitted  symbols.  For  BPSK  symbols, 
the  detection  is  reduced  to 


bk  =  sign(Re{w^0yfc}), 

while  for  4-QAM  symbols,  the  detection  can  be  elaborated  as 
h  =  sign  (Re  {w^Gyfe})  +  sign  (Im  {w^cyfe})  j. 


in  which  j  =  ^/— I.  The  combining  weight  vector  wk,0  in  (10) 
of  size  2 L  is  chosen  to  maximize  the  SINR  at  the  destination 
Dk  which  is  given  by 


SIMR(w^)  = 


EUw»  ,fcC  ibi  + 


(11) 


i.e., 


w ko  —  arg  maxSINR(wju). 
w  k 

Note  that  in  (1 1),  the  expectation  in  the  numerator  is  taken  over 
the  random  variable  bk,  while  the  expectation  in  the  denomi¬ 
nator  is  taken  over  the  random  variables  bi ,  l  ^  k  and  all  inde¬ 
pendent  noise  terms  in  . 

III.  System  Performance  Analysis 

To  determine  the  maximum  SINR  weight  vector  wk,0  and  the 
corresponding  SINR  at  the  destination  Dk,  we  first  define 


K 

u *=  E  Hac*cfHffc  +  rfc, 

1=1, l^k 


Tk±E{nkn*} 


(IL  0  \ 

Vo  Pra2Ur,dkU^dk+ILJ- 


At  each  destination  Dk ,  it  combines  the  signals  received  in 
Phase  1  and  the  signals  received  from  the  relay  in  Phase  2  to 
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Then  the  SINR  in  (11)  with  any  given  combining  weight  vector  and  denote  by  C'-1  the  fcth  column  vector  of  the  matrix  . 
Wfc  can  be  written  as  Then,  after  some  algebraic  calculations,  we  can  see  that 


SINR(wfc)  = 


wfHfc,fccfc| 

wfUfcWfe 


(12) 


It  is  easy  to  check  that  U/,;  is  Hermitian,  so  are  u[l,/2)  and 
U-(1/2)d  According  to  Schwartz  inequality,  we  have 


wfHfe;fcc  k 

wfUfeWfe 


< 


wfe  U|Ufe 


wfuIuIwfe 


wfuI 


u.  2H k,k£k 


wfui 


where  the  equality  holds  when  w^U^1//2^  = 

Thus,  the  maximum  SINR  weight 
vector  wjfejG  is  given  by 


w  u 


=  U.  1Hfc  fcCfc. 


The  corresponding  maximum  SINR  at  the  destination  Dk  with 
the  optimum  weight  vector  Wfc?G  is  equal  to 

SINRfc  =  cf  H£feUfc-lHMcfe.  (13) 

In  order  to  optimally  allocate  power  to  all  sources, 
we  further  exploit  SINRfc  in  (13)  as  follows.  We  de- 

fine  Cfe  =  diag(ci,...,cfe_i,cfc+i,...,cic),  which  is  a 
L(K  -  1)  x  (K  —  1)  block  diagonal  matrix  formed  by 
placing  all  code  vectors  except  Cfc  in  the  diagonal  positions, 

and  Hfc  =  . . . ,  Hfe+i,fe, . . . ,  H K}k\  which  is  a 

2  L  x  L(K  —  1)  interference  channel  matrix.  Using  the  above 
notation,  SINRfc  can  be  expressed  as 


SINRfc  =  cf  (HfcCfcCf  Hf  +  rfe)  *  Hfc.fcCfc. 

(14) 

According  to  the  Woodbury  matrix  inversion  lemma  [27],  we 
have  (HfcCfcCf  Hf  +  Tfc)-1  =  r* 1  -  ^Cfc^-i  + 
CfHfr^^fcCfc)-1^^1. 

Let  us  now  define  the  following  matrix  F^)  =  (  f where 


/ife)  =  VP^Pn  [c*HfmidkHSn4kcn 

+  (Pra2Ur4kU^dk  +  IL)  1 

^  Hr?c^ HSri?rcn]  (15) 

'We  may  represent  U  k  in  terms  of  its  eigenvalues  A, ,  A2 . XK  and  their 

corresponding  eigenvectors  q , .  q2 . q  a  as:  U k  —  -Vqzq/7,  where 

A/  >  0  for  1  <  /  <  Ii  and  qfqz/  .=  ^  V .  Then,  U^1/2)  = 

Az(1/2)q;q/H  and  u7(1/2)  =  V-(1/2)qzqf . Xt  is  easy t0  check 

that  both  U^1/2)  and  are  Hermitian. 


(k) 

where  ft  y  contains  the  channels  and  the  cross-correlation  be- 

k 

tween  the  intended  source  Sk  and  the  interfering  sources,  while 
contains  the  channels  and  the  cross-correlation  among  in¬ 
terfering  sources.  Based  on  (14)— (16),  we  can  represent  SINRfc 
as 


SINRfc  =  fif  -ff)H(lK_1+Ff)  (17) 


where  the  superscript  ( •  indicates  the  corresponding  desti¬ 
nation  Dk . 

Note  that  if  the  channels  are  quasi-static  (i.e.,  constant) 
over  each  information  symbol  period,  i.e.,  dk  =  hSn,dk, 
hillr  =  hSn^  and  h^d%  =  hr4k  Vi  =  1,2,  .  and 
Vn,  k  =  1,  2, . . . ,  K,  then  in  (15)  is  reduced  to 


frnn  ~  Pmn\f^mPn 


x  Km,dkh^4k  + 


ol  Pr  \hr  dk  | 


a2Pr\h 


r,dk  | 


!  +  l 


(18) 


As  an  illustration  example,  when  there  are  only  two  source-des¬ 
tination  pairs  with  one  relay  in  the  network,  i.e.,  AT  =  2,  the 
SINR/,  in  (17)  with  quasi-static  channels  can  be  specified  as 


SINRfc  =  fj» 


Ak) 
J  kj 


(i  +  Zf)"1 


—  Pk  \hSk,dk\  + 


2  a2Pr  \hr,dk\2  \hs 


Oi2Pr  \hr,dk\  +  1 


p2PkPj 


^sk  Ak 


2Pr\hr^.\Zh:  h„ 


a2Pr\hr,dk\2  +  l 


1  +  Pj  \hSj^dk  |  + 


<x2Pr\hr,dk\‘'Pj\haj,r 

^Pr\hr,dJ2+l 


where  p  —  p\2  —  P21  and  the  subscript  j  =  2  if  k  =  1  while 
j  =  1  if  k  =  2.  Furthermore,  when  the  code/channel  vectors 
Cfc  are  orthogonal  to  each  other,  i.e.,  pkj  =  c[ Cj  =  0,  for 
j  /  k ,  then  =  0  for  quasi-static  channels.  Hence  the  SINR 

at  destination  Dk  is  given  by 


SINRfc  = 


n  It  |2  0?Pr  \hr  dk  |  Pfc  |Zlss.  r| 

=  Pk\riSk,dk\  + 


a2Pr\hr,dk\  +1 


(19) 


which  is  the  sum  of  the  SINRs  of  the  direct  link  and  the  relay 
link. 


IV.  Optimum  Power  Assignment  Under  SINR 
Constraints  for  All  Source-Destination  Pairs 

In  this  section,  we  determine  the  optimum  power  assignment 
for  the  sources  and  the  relay  that  minimizes  the  total  power  con- 
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sumption  under  the  condition  that  the  SINR  requirement  of  each 
source-destination  pair  is  satisfied.  First,  we  consider  the  power 
optimization  for  the  relay  network  in  a  general  setting  where 
codes/channels  of  different  sources  may  have  arbitrary  correla¬ 
tion.  Then,  we  discuss  a  simplified  power  optimization  scheme 
for  a  special  case  where  codes/channels  of  different  sources  are 
orthogonal,  and  provide  an  intuitive  interpretation  for  the  pro¬ 
posed  scheme. 

A.  Minimization  of  Total  Power  Consumption 

Let  us  assume  that  the  SINR  requirement  for  the  source-desti¬ 
nation  pair  (5fc,  Dk)  is  7*.,  fc  =  1, 2, . . . ,  K.  Then,  the  problem 
of  optimizing  power  to  minimize  the  total  power  consumption 
and  satisfy  all  source-destination  SINR  requirements  can  be  for¬ 
mulated  as 


—  det  ( Iff-i  +  FE 


(1  +  SINRfc) 


where  the  last  equality  follows  from  the  expression  of  SINR/. 
in  (17).  Thus,  we  have 


det(lK  +  F(fe)) 
det  (!*_!+ F^' 


1  +  SINRfe  = 


We  note  that  for  moderate  or  high  SINR,  we  may  approximate 
1  +  SINRfc  «  SINRfc  and  1  +  «  f^k\  Vfc,  l  =  1, . . . ,  K. 

So,  from  (26)  we  may  approximate  SINR/,  as  follows: 


SINRfc 


det  (F(fe))  _  P/tdet  (G(fe)) 


det  F 


det  G 


/  min P1,„.,PK;Pr  Efe= l  A  +  P  (20) 

l  s.t.  SINRfe  >  7fe,  1  <  k  <  K 

where  the  transmission  power  terms  Pi ,  P2 , . . . ,  Pk  and  Pr  are 
all  nonnegative. 

Let  us  define  an  auxiliary  parameter,  which  can  be  viewed  as 
a  normalized  power  factor  at  the  relay,  as  follows: 

W  =  a2Pr  (21) 

where  a  is  the  amplification  factor  specified  in  (8).  The  auxiliary 
parameter  x  will  play  a  key  role  in  the  optimization  procedure. 
Let  us  also  denote  by  =  (gmn)  a  matrix  with  elements 

(k)  A  -HttH  tt  _i_  HttH  tt// 

i/mn  -n-'Sri  "r  ,r*-*-r.dk 

x  (xUrAU^dk+IL)  1  Hr:dkUSn:rcn  (22) 

for  any  m,  n  =  1,  2, . . . ,  iT.  Then,  from  (15),  we  can  represent 
each  entry  in  by 

f(k )  —  A/p  P  n(k) 

Jmn  ~  V  r mr ny-mri' 

It  is  straightforward  to  verify  that 


det  (fW)  =  (  npM  det  (G(fc)) 


det  F 


Ak)\  — 


1=1, i^k 


II  **K 


From  (17),  we  know  that  Ik-i  +  F^;  is  in 
on  the  Schur  complement  formula2,  we  have 


det  (lK  +  F<fc>) 


1  -I- 

1  +  Jkk 


=  det  (Ik-i 


1  1  f(fc)  f  (k)H 

L^hk  h 


2If  matrix  D  is  invertible,  then  [27] 


Ik-i  +  Ff 


Ik-i+F^ 


=  det(D)  •  det(A  -  BD_1C) 


The  above  approximation  is  asymptotically  tight  for  high  SINR. 
Based  on  this  approximation,  the  optimization  problem  in  (20) 
can  be  written  as 


mm  P1,...,PKiPr 


ELi  Pk  +  P, 

Pk  det(G<fc>) 
det.(Gifc)l  _ 


>  Ik,  1  <k<K. 


F(P  is  invertible.  So,  based 

k 


Let  V  denote  the  set  of  feasible  solutions  for  the  optimization 
problem  in  (28),  i.e. 

V  =  {Pi,...,Pjc,Pr  ISINRfc  >  7fe,  VI  <  k  <  K}. 

We  may  further  partition  the  set  V  into  disjoint  subsets  such  as 
V=  \Jvx 


where 

V  (p  p  p  SINRfe  >  7fc,  VI  <  k  <  K, 
Vx  =  iPp...,PK,Pr  a2Pr  =  x 


for  any  x  >  0. 

We  note  that  for  any  given  value  of  the  auxiliary  parameter  x 
in  (21),  according  to  (8)  the  transmission  power  at  the  relay  Pr 
can  be  determined  as 

K 

Pr=  A  =XP  PkPsk,r  +  XL.  (29) 

k= l 

Thus,  for  any  given  x  >  0,  the  optimization  problem  in  (28) 
over  the  feasible  set  V,  becomes 


min/> . ph  Efc=i  +  1)  Pfe  + 


jkdet(G^] 

det(G<») 


,  1  <  k  <  K. 


We  observe  that  in  (30),  for  any  given  x  >  0, 

(7/,det(G^))/(det(G^))  is  a  constant  which  is  in¬ 
dependent  of  Pk  (k  =  1,  2, . . . ,  K).  Hence  the  minimal  value 
in  (30)  is  obtained  when  all  constraints  hold  with  equality,  i.e. 

7fcdet  (g^) 

P*=  det  (GW)  '  (31) 
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Then,  the  corresponding  minimal  total  power  of  (30)  is 

A  K  det  (Gf)) 

v (x)  =  Lj  ^  det  (GW)  +  XL  (32) 

which  is  a  function  of  x  >  0.  Let  us  denote  v*  as  the  minimal 
value  of  the  objective  function  over  (28)  in  the  feasible  set  V. 
Then,  we  can  see  that 


v*  =  min  !;(#).  (33) 

tc>0 

The  above  discussion  shows  that  we  are  able  to  convert  the 
optimization  problem  in  (20)  over  a  multidimension  space  to 
the  minimization  problem  in  (33),  which  depends  only  on  one 
variable  x  >  0,  i.e.,  over  a  one-dimension  space.  The  minimiza¬ 
tion  of  v{x)  in  (32)  can  be  easily  solved  by  a  numerical  search 
for  the  optimal  value  of  the  parameter  x  >  0.  With  the  optimal 
value  x *  that  minimizes  the  function  v(x)  in  (32),  we  can  obtain 
the  corresponding  optimal  power  and  P*  based  on  (31)  and 
(29),  respectively. 

B.  Simplified  Optimization  With  Orthogonal  Codes/Channels 

In  the  previous  subsection,  we  solved  the  optimization 
problem  for  a  general  multisource  multi-destination  relay 
network  with  arbitrary  correlation  among  user  codes.  In  this 
subsection,  we  are  able  to  further  simplify  the  optimization 
when  the  signatures  of  different  sources  are  orthogonal  and  the 
fading  channels  are  quasi-static  during  the  transmission  period 
of  a  signature  code. 

In  particular,  for  multi- source  multi-destination  relay  net¬ 
works  with  orthogonal  transmissions,  the  cross-correlation 
matrix  R  in  (1)  is  an  identity  matrix,  so  both  and  G ^  in 
(27)  are  diagonal  matrices,  and 

det  K  0  =  njWjL  =  (34) 

n£ia?! 


where  gj^  is  specified  in  (22).  With  the  assumption  of  quasi¬ 
static  channels,  we  can  write  h$,r  =  hahft,  h^d  =  h3kidk, 
and  h^\k  =  hr^k ,  Vz  =  1, 2, . . . ,  L,  and  Vfc  =  1, 2, . . . ,  K. 
We  note  that  in  the  minimal  total  power  in  (32),  fi3k,r  =  ^  r, 
where  (j2Sk  r  is  the  variance  of  the  source-relay  channel  hSk^r, 
i.e.,  (Jgk  r  =  E{\hSkir\2}.  Thus,  by  substituting  (34)  into  (32), 
we  have 


■j(X)  =  jr  lk  ^ r  + 


*=1  KaI2  +  ^ ^ 


+  xL 


K 


k= 1 


=  7  k[ck  +  hx  + 


ak 


x  +  dk 


(35) 


where 


fi'k  — 


I  hsk,r 


A2\hr4k\ 


1  - 


dk  I 


|2  2 


A  \  hr4k  | 


bk  = 


+ 


Klk  A 
_  1  |  hSk. 

A  A2\hr4k\ 


2  a2 

sk,r 


\hsm 


dk  I 


dk  = 

A\K,dkr 

and  A  =  \h3kjdk  |2  +  |^sfc,r|2- 


(36) 


We  note  that  for  any  k  =  1,2,...,  K,  if  a/,;  >  0,  each  term 
Ck  +  bkX-\-(ak)/(x-\-dk)  in  (35)  is  convex  with  respect  to  x  >  0, 
and  it  can  be  minimized  by 


x  =  max 


(o,  —dk  + 


If  ak  <  0 ,  Cfc  +  bkX  +  (a/,.)/ (x  +  dk)  is  increasing  in  [0,  +oo),  it 
implies  that  the  minimum  point  occurs  at  x  =  0.  Let  us  denote 

^min  =  mm(x1,x2,...,xK)  (37) 

^max  =  ma,x(xi1x2,..-,xK)  (38) 


where 


max  (o,  -dk  +  yj §£)  if  ak  >  0;  ^ 

0  if  ak  <  0; 


for  any  k  —  1,  2, . . . ,  K.  Then,  each  term  Ck  +  bkX  +  (ak)/ (x  + 
dk)  in  (35)  is  decreasing  in  the  range  [0,  xmin]  and  increasing  in 
the  range  [rrmax,  +oo),  which  implies  that  the  optimal  solution 
x *  that  minimizes  the  function  v(x)  in  (35)  is  bounded  as 


<  X*  <  Xr, 


(40) 


Thus,  to  find  the  optimal  solution  x *,  we  only  need  to  search 
within  the  range  [xm[ni  rrmax].  We  note  that  xm[n  =  0  if  there 
exists  k  £  {1,  2, ... ,  K}  such  that  ak  <  0. 

From  (39),  we  can  see  that  the  necessary  condition  for  Xk  >  0 
is 


>  bkd2k , 


which  implies  that 


l^sjyrl  Vlr,dk  I 


C 7% 


>  |  h 


s k  | 


(41) 


In  a  nonfading  or  slow -fading  scenario,  where  the  coherence 
time  of  the  channel  is  much  longer  than  the  delay  requirement 
of  the  channel,  the  channel  gain  remains  roughly  constant  over  a 
period  of  processing  time  (for  example,  transmitting  several  data 
packets)  [28].  In  this  case,  we  may  safely  assume  that  \h{j  |2  « 
af  j  during  the  processing  period,  where  afj  is  the  variance 
of  the  channel  hij.  With  such  an  approximation,  the  necessary 
condition  in  (41)  is  reduced  to 


<Jr,dk  >  ask,dk  *  ^2) 

If  the  necessary  condition  in  (42)  does  not  hold,  i.e.,  a^dk  < 
< 72Sk  dk  for  all  k  =  1,2,...,  K ,  then  xmax  =  0,  which  implies 
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that  the  optimal  solution  x*  =  0.  So,  the  corresponding  op¬ 
timum  power  at  the  relay  is  P*  =  0,  which  means  that  the  relay 
is  not  needed  in  this  case.  In  other  words,  if  each  relay-destina¬ 
tion  channel  link  is  weaker  than  the  intended  source-destination 
channel  link,  then  the  use  of  the  relay  is  not  helpful. 

In  addition,  if  ak  <  0,  Vfc  =  1,  2, . . . ,  K,  then  each  term 
Ck  +  bkX  +  (ak)/(x  +  dk)  in  (35)  is  an  increasing  function  of 
x  E  [0,  -|-oo),  therefore  the  optimal  solution  is  x*  =  0.  Given 
that  ak  <  0  is  equivalent  to 

\h'Sk-Ak  I  ask,r  —  (\hm,dk\  +  )  \hr,dk  |  (43) 

it  implies  that  if  all  channels  \hSk^k  |2  (k  =  1,  2, . . . ,  K)  that 
correspond  to  direct  links  are  strong  enough  or  all  channels 
of  the  relay  links  \hr^k  |2  {k  =  1, 2, . . . ,  K)  are  weak,  then 
avoiding  use  of  the  relay  can  leads  to  power  savings.  In  other 
words,  if  there  is  at  least  one  relay-destination  link  that  is  better 
than  the  corresponding  source-destination  link,  the  relay  should 
forward  its  received  signal  to  the  destinations  in  order  to  help 
reduce  the  overall  power  consumption  while  maintaining  at  the 
same  time  the  target  performance  level. 

Finally,  for  a  symmetric  system  where  g2Si  di  =  a2  do  = 

•••  =  =  aLr  =  •••  =  AK,r  and  "r’d,  = 

(T2d2  —  •  •  •  =  (r2d ,  we  are  able  to  obtain  an  analytical  solu¬ 
tion  for  the  minimization  of  v(x).  In  this  case,  all  dk  are  equal, 
denoted  by  d  =  dk,  Vfc.  Then 

v(x)  =  c+bx-\ - (44) 

x  +  d 

where  c  =  Efc=i7fcCfc,&  =  i'7kfa,a  =  Efc=i7fcafc- 

Consequently,  the  value  x *  that  minimizes  v(x)  in  (44)  is 

—  I  ~d  +  if  a  >  0,  (45^ 

\0  if  a  <  0. 

Substituting  (45)  into  (31)  and  (29),  we  obtain  the  optimal  power 
PjJ  and  P*,  respectively. 

We  note  that  for  a  symmetric  system,  all  ak  are  equal.  When 
the  SINR  requirement  is  high  enough  at  each  destination, 
L/(K^k)  is  relatively  small  compared  to  (a2k  r)/(A),  so  all 
bk  s  are  approximately  the  same  (h  =  l,2,...,iT).Asa  con¬ 
sequence,  the  ratio  a/ bin  (45)  is  independent  of  7*..  Hence,  the 
optimal  value  x*  in  (45)  is  independent  of  the  SINR  thresholds 
7fc  in  this  case. 


C.  Equal  Power  Allocation  Scheme 

In  this  subsection,  for  comparison  purposes,  we  discuss  the 
equal  power  allocation  scheme  where  all  the  sources  and  the 
relay  are  allocated  the  same  power  P.  In  this  case,  the  corre¬ 
sponding  parameter  x  in  (21)  is 


x  = 


1 


E 


K 

k= 1 


a 


2 

skX 


(46) 


To  find  the  optimal  power  P  that  minimizes  the  total  power 
of  the  system  under  the  constraints  that  the  SINR  requirements 
of  all  source-destination  pairs  are  satisfied,  we  can  follow  the 
procedure  described  in  the  previous  subsection,  i.e. 


minp  PEfc=iC 


xat 


s .  t . 


7fcdet(G<-fc)) 
-  det(G(fc)) 


+  l)  +  xL, 

,  k  =  1,2, 


,  K. 


(47) 


We  note  that  the  objective  function  in  (47)  is  increasing  in  terms 
of  increasing  P  >  0.  Thus,  the  optimal  solution  P*  of  (47)  is 
given  as 


p* 


max 

k=l,...,K 


7fcdet  (Gg0)  I 
det  (G(fc))  J 


(48) 


Since  the  channel  quality  is  not  the  same,  in  general,  for  all  links, 
it  is  implied  that  the  terms  (7/cdet(G^))/(det(G^^))),  k  = 
1.2,...,  K  in  (47)  are  not  equal  in  general,  which  means  that 
the  equality  in  (47)  might  not  hold  for  all  L  As  a  result,  the 
equal  power  allocation  strategy  generally  spends  more  power 
than  what  is  needed. 


V.  Optimum  Power  Assignment  Under  a  Total  Power 
Budget  Constraint 

In  this  section,  we  apply  the  methodology  that  we  developed 
in  Section  IV  to  solve  a  related  max-min  SINR  based  power 
optimization  problem  for  multi- source  multi-destination  relay 
networks.  Specifically,  we  design  a  power  assignment  scheme 
that  maximizes  the  minimum  SINR  among  all  source-destina¬ 
tion  pairs  subject  to  any  given  total  power  budget.  We  note  that 
such  a  scheme  introduces  a  type  of  fairness  among  different 
source-destination  pairs.  The  optimization  problem  can  be  for¬ 
mulated  as  follows: 

(  maxPl . pK;Pr  minfc=ii...!K{SINRfe} 

t  s  •  t  •  Ef=l  A  +  Pr  <  Aotal 

where  Ptotai  is  the  given  power  budget  of  the  system. 

Let  us  denote 

x  =  minfc=ii...iK-{SINRfc} 

then  (49)  is  equivalent  to 

ma  xpli...ipA.;pr  z, 

S.t.  Efc='i  A  •  Pr  5~  -Ptotah 

SINR/.  >#fk  =  l,...,K. 

According  to  (27),  we  may  approximate  SINR*,  as 


(49) 


(50) 


(51) 


SINRfc 


Pk det  (G(fe)) 
det  (Gf) 


(52) 


where  the  matrices  G  W  and  G^  are  specified  in  (22).  If  we 

define  an  auxiliary  parameter  x  =  a2Pr ,  then  for  any  given 
x  >  0,  the  transmission  power  at  the  relay  Pr  is  determined  as 

K 

Pr  =  =  X  Pkftsk  ,r  +  XL.  (53) 

or  f — ^ 


As  a  consequence,  we  can  reformulate  the  problem  in  (51)  as 

(  maxplv..?px  My  ^ 

■  +  1)  Pk,  +  xL  <  Ptotai 

i«(gT) 


S.t. 


Et;=l  (X0Sk-,r 


Pk>z- 


det(G(fc))  1 


1  <  /,:  <  K. 


(54) 


7 


CHEN  et  al.\  JOINT  POWER  OPTIMIZATION 


2377 


For  any  x  >  0,  we  denote  the  maximal  value  of  the  objective 
function  in  (54)  as  z(x).  Then,  z(x)  must  satisfy 

det  (g^) 

Pk  >  z(x)  ■  —  >  Vk  =  l,2,...,K.  (55) 

Substituting  the  above  constraints  into  the  total  power  constraint 
in  (54),  we  have 

k  det  (G^) 

53  (x&tnr  +  1)  z(x)  •  det^G(fc)^  +xL<  Ptotal.  (56) 


Thus,  for  any  fixed  x  >  0,  the  maximal  value  of  the  objective 
function  in  (54)  is 


z(x)  = 


Ptotal  xL 


£ 


K 

fe= 1 


{xBSk,r 


+  !)• 


det(G(fc)) 


(57) 


We  can  see  that  the  maximization  of  z(x)  can  be  implemented 
by  a  numerical  search  of  single  (one-dimension)  parameter  (x  > 
0).  If  we  denote  the  optimal  parameter  of  (57)  as  x *,  then  the 
optimum  power  assignment  to  the  sources  is 


optimization  problem  in  (49)  with  an  equal  power  assignment 
can  be  expressed  as 


maxp 

S.t.  (K  +  1  )P  <  P total 5 


(60) 


where  all  the  sources  and  the  relay  are  assigned  the  same  trans¬ 
mission  power  P.  With  the  SINR  approximation  in  (52)  and 
M  =  mini<fc<ir{SIMRfc},  we  have 


maxp 

s.t. 


P  < 


K+ 1 


P>Z  • 


det(G<C) 

det(G(fc))  ’ 


k  =  l,..^K. 


(61) 


Thus,  the  optimal  solution  is  P  —  (Ptotal) /(^  +  1)  and  the 
maximal  value  of  the  worst-case  SINR  is  given  by 


z  — 


Ptotal 

•  maxi <k<K 


det  (g)  ) 

det  (g<*>)  j  ' 


(62) 


We  can  see  that  in  the  equal  power  assignment  scheme,  the  min¬ 
imum  SINR  value  among  all  source-destination  pairs  is  directly 
determined  by  the  given  total  power  budget  value. 


det  (g^) 

Pt‘  =  *(X1 '  det  (GW)  ’  k  =  1’2--K  <58) 

and  the  corresponding  power  assigned  to  the  relay  is 

P*  =  x*  (jlPk^r  +  Lj.  (59) 

We  note  that  z(x)  in  (57)  should  be  nonnegative  for  any  given 
x,  which  implies  that  Ptotal  ~  xL  >  0,  i.e.,  x  <  (Ptotai)/(P). 
Thus,  we  only  need  to  search  the  single  variable  x  over  the  in¬ 
terval  [0,  (Ptotai)/(£)]  t0  obtain  the  optimal  solution  x*  that 
maximizes  the  function  z(x). 

From  (52)  and  (58),  we  can  see  that  for  any  k  —  1,  2, . . . ,  K, 
SINRfc  =  z(x*),  which  implies  that  the  optimum  power  as¬ 
signment  achieved  through  the  max-min  SINR  based  optimiza¬ 
tion  in  (49)  leads  to  the  same  SINR  values  for  all  source-des¬ 
tination  pairs.  We  recall  that  for  the  total  power  minimization 
problem  in  (20),  the  optimal  power  assignment  ensures  that  the 
resulting  SINR  for  each  source-destination  pair  is  equal  to  the 
SINR  requirement  of  the  corresponding  source-destination  pair, 
i.e.,  SINR*;  =  7*.,  k  =  1, . . . ,  K.  It  is  natural  to  expect  that 
when  all  source-destination  SINR  requirements  7*.  in  (20)  are 
the  same,  the  power  assignment  strategy  for  the  max-min  SINR 
based  optimization  and  that  for  the  total  power  minimization 
based  optimization  should  be  the  same.  This  intuitive  interpre¬ 
tation  can  be  seen  from  the  derivations  in  (32)  and  (57).  We  note 
that  for  reasonably  high  power  budget  Ptotai,  the  term  xL  in  (57) 
can  be  ignored,  and  the  maximization  of  z(x)  in  (57)  is  equiv¬ 
alent  to  the  minimization  of  the  dominator  which  is  related  to 
the  objective  function  in  (32).  If  all  the  source-destination  pairs’ 
SINR  requirements  in  (32)  are  the  same,  i.e.,  7*.,  k  =  1 .....  /T, 
then  the  optimal  solutions  x*  of  z(x)  in  (57)  is  also  optimal  in 
minimizing  v(x)  in  (32). 

For  comparison  purposes,  in  the  following  we  illustrate  the 
max-min  SINR  problem  under  an  equal  power  allocation.  The 


VI.  Numerical  Results 

In  this  section,  we  perform  numerical  studies  to  illustrate  the 
proposed  optimum  power  assignment  algorithms.  In  our  studies, 
we  consider  a  slow  fading  scenario  where  the  coherence  time  of 
the  channel  is  much  longer  than  the  delay  requirement  of  the 
channel  and  then  approximate  \hij\2  by  a2  -  [28].  The  channel 
gain  for  each  channel  link  is  assumed  to  follow  a  path  loss 
model,  where  the  variance  of  channel  coefficient  is  given  by 
a2  j  =  Sfj  (fj  G  {sfc,  dk,  r})  with  6ij  as  the  distance  of  the 
channel  link  and  A  as  the  path-loss  exponent  (A  =  3  in  our  nu¬ 
merical  studies). 

A.  Total  Power  Minimization 

In  the  first  set  of  numerical  studies,  we  illustrate  the  optimum 
power  assignment  that  minimizes  the  total  power  consumption 
under  the  condition  that  the  SINR  requirement  of  each  source- 
destination  pair  is  satisfied.  First,  we  consider  a  system  with 
two  source-destination  pairs  and  one  relay,  i.e.,  K  =  2,  and  the 
cross-correlation  of  the  two  source  codes  is  p  =  0.25.  We  study 
two  cases  of  SINR  requirements  for  the  two  source-destination 
pairs:  (i)  [71,72]  =  [10, 10]  dB,  and  (ii)  [71,72]  =  [10,20]  dB. 

Fig.  2  plots  the  total  power  consumption  with  varying  param¬ 
eter  x  >  0  for  an  asymmetric  system,  where  the  values  of  the 
distance  between  nodes  are  set  as  follows:  =  2, 8Sl^  — 

=  =  8 s 2 , <7 1  =  =  1}  ^.82,^2  =  ^  <ffld 

Sr  d2  =  2.  The  optimal  values  of  the  parameter  that  minimize 
the  total  power  consumption  are  x*  =  0.63  and  x*  =  1.52 
for  [71,72]  =  [10, 10]  dB  and  [10,20]  dB,  respectively.  Based 
on  the  optimal  value  x *,  we  obtain  the  corresponding  optimal 
power  assignment  Pi  ^2  and  Pr  according  to  (31)  and  (29), 
as  listed  in  Table  I.  In  this  table,  we  also  compare  the  optimal 
power  values  obtained  by  our  proposed  approximation  method 
and  those  obtained  by  exhaustive  search  based  on  the  optimiza¬ 
tion  in  (20).  We  observe  that  the  optimal  power  values  obtained 
by  the  two  methods  are  almost  indistinguishable.  In  Fig.  3,  we 
plot  the  total  power  consumption  with  varying  parameter  x  >  0 
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Total  power  vs.  variable  x,  K=2 


Fig.  2.  Minimization  of  the  total  power  consumption  under  varying  parameter 
x  for  an  asymmetric  multi-source  multi-destination  relay  network  with  given 
SINR  constraints  (K  —  2). 


Total  power  vs.  variable  x,  K=2 


Fig.  3.  Minimization  of  the  total  power  consumption  with  varying  parameter  x 
for  a  symmetric  multi-source  multi-destination  relay  network  under  given  SINR 
constraints  ( K  —  2). 


for  a  symmetric  system,  where  the  values  of  the  distance  be¬ 
tween  nodes  are  set  as:  8Sudl  =  8S2,d2  =  8Sltd2  =  8S2,dl  =  2 
and  <5Sl?r  =  <5S2?r  =  8r^dl  =  8r^d2  =  1.  The  optimal  values 
of  the  parameter  are  x*  =  0.32  and  x*  =  0.59  for  [71,72]  = 
[10, 10]  dB  and  [10,  20]  dB,  respectively.  Based  on  the  optimal 
value  x *,  we  obtain  the  corresponding  optimal  power  alloca¬ 
tion  Pi,P2  and  Pr  again  according  to  (31)  and  (29),  as  listed 
in  Table  II.  We  observe  that  when  the  SINR  requirements  of  the 
two  source-destination  pairs  are  equal,  the  power  assignments  to 
the  two  sources  are  the  same.  However,  when  the  SINR  require¬ 
ments  are  different  ([71,72]  =  [10,20]  dB),  the  equal  power 
assignment  is  not  optimum  anymore.  We  can  see  in  Table  II  that 
the  optimal  power  values  obtained  by  our  proposed  approxima¬ 
tion  method  are  indistinguishable  from  those  obtained  by  the 
exhaustive  search  based  on  the  optimization  in  (20). 


Total  power  vs.  variable  x,  K=3 


Fig.  4.  Minimization  of  the  total  power  consumption  with  varying  parameter 
x  for  an  asymmetric  multi-source  multi-destination  relay  network  under  given 
SINR  constraints  (K  —  3  ). 


Total  power  vs.  variable  x,  K=3 


Fig.  5.  Minimization  of  the  total  power  consumption  with  varying  parameter  x 
for  a  symmetric  multi-source  multi-destination  relay  network  under  given  SINR 
constraints  ( K  —  3). 


We  now  repeat  our  studies  for  a  system  with  three  source-des¬ 
tination  pairs  and  one  relay,  i.e.,  K  =  3,  where  the  signa¬ 
tures  of  the  three  sources  are  orthogonal  to  each  other.  We  ex¬ 
amine  both  an  asymmetric  case  (Fig.  4)  and  a  symmetric  case 
(Fig.  5).  In  both  figures,  we  consider  three  sets  of  SINR  re¬ 
quirements  for  the  three  source-destination  pairs,  namely:  (i) 
[71,72,73]  =  [20,20,20]  dB,  (ii)  [71,72,73]  =  [10,20,20] 
dB,  and  (iii)  [71,72,73]  =  [10, 10,20]  dB.  In  Fig.  4,  we  con¬ 
sider  an  asymmetric  system  where  the  distance  values  are  set  as 
$sk,r  =  1?  (&  =  I?  2, 3),  8r,dl  =  1, 8rdi2  =  2, 8r^  =  3  and 
8Sudl  =  2, 8S2^d2  =  3  and  8S 3^ds  =  4.  In  this  case,  the  optimal 
values  of  the  parameter  that  minimize  the  total  power  consump¬ 
tion  are  x*  =  2.34,  x*  =  2.77  and  x*  =  3.17  for  the  three 
sets  of  SINR  requirements,  respectively.  Based  on  the  optimal 
values,  the  optimum  power  assignments  Pi,  P2,  P3  and  Pr  can 
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TABLE  I 

Optimum  Power  Assignment  That  Minimizes  the  Total  Power  Under 
Given  SINR  Constraints  Using  the  Proposed  Method  and  the 
Exhaustive  Search  Method  for  an  Asymmetric  Setting  (K  =  2) 


(7L72)(dB) 

Pi(dB) 

F2(dB) 

Pr(dB) 

Proposed  method 

(10,10) 

13.14 

19.89 

20.54 

Exhaustive  search 

(10,10) 

13.22 

19.91 

20.41 

Proposed  method 

(10,20) 

11.59 

27.34 

29.74 

Exhaustive  search 

(10,20) 

11.76 

27.28 

29.70 

TABLE  II 

Optimum  Power  Assignment  that  Minimizes  the  Total  Power  Under 
Given  SINR  Constraints  Using  the  Proposed  Method  and  the 
Exhaustive  Search  Method  for  a  Symmetric  Setting  (K  =  2) 


(7i,72)(dB) 

Pi(dB) 

P2(dB) 

Pr(dB) 

Proposed  method 

(10,10) 

14.56 

14.56 

15.86 

Exhaustive  search 

(10,10) 

14.91 

14.91 

15.31 

Proposed  method 

(10,20) 

13.25 

23.25 

22.44 

Exhaustive  search 

(10,20) 

13.32 

23.32 

22.34 

TABLE  III 

Comparison  of  the  Total  Power  Consumption  Resulting  from  the 
Optimum  Power  Assignment  Scheme  and  the  Equal  Power  Assignment 
Scheme  (K  =  3) 


71  >  72 , 73  (dB) 

Total  Power  Consumption  (dB) 

Optimum 

Equal 

Asymmetric 

System 

[20,20,20] 

37.15 

41.58 

[10,20,20] 

36.82 

41.58 

[10,10,20] 

35.63 

41.58 

Symmetric 

System 

[20,20,20] 

29.79 

30.30 

[10,20,20] 

28.25 

30.30 

[10,10,20] 

25.85 

30.30 

be  obtained  accordingly  based  on  (31)  and  (29).  In  Fig.  5,  we 
consider  a  symmetric  system  with  SSk,r  =  Sr,dh  =  1  ,(fc  = 
1,2,3)  and  SSk,dk  =  2,  (k  =  1,2,3).  The  optimal  values  of 
the  parameter  are  x*  =  0.77,  x*  =  0.76  and  x*  =  0.75  for 
the  three  sets  of  SINR  requirements,  respectively.  We  can  see 
that  the  optimal  values  x*  in  the  symmetric  system  are  almost 
the  same  for  the  three  sets  of  SINR  requirements.  This  result  is 
consistent  with  the  theoretical  discussion  that  the  optimal  value 
x *  is  independent  of  the  SINR  requirement  7^  but  it  depends  on 
the  channel  condition. 

In  Table  III,  we  show  the  power  efficiency  of  the  proposed 
optimum  power  assignment  scheme  by  listing  the  resulting  total 
power  consumption,  compared  to  the  total  power  consumption 
that  results  from  the  equal  power  assignment  scheme.  We  con¬ 
sider  both  an  asymmetric  case  (the  system  setup  is  the  same  as 
that  in  Fig.  4)  and  a  symmetric  case  (the  system  setup  is  the  same 
as  that  in  Fig.  5).  We  can  see  that  in  the  asymmetric  system, 
the  power  savings  of  the  optimum  power  assignment  scheme  is 
4-6  dB  for  the  three  sets  of  SINR  requirements.  The  more  un¬ 
balanced  the  SINR  requirements  of  the  three  source-destination 
pairs  are,  the  more  power  savings  of  the  optimum  power  assign¬ 
ment  scheme  compared  to  the  equal  power  assignment  scheme 
can  be  achieved.  In  the  symmetric  system,  the  power  savings 
of  the  optimum  power  assignment  scheme  is  0.5-4.5dB  for  the 
three  sets  of  SINR  requirements.  Comparing  the  results  between 
the  asymmetric  system  and  the  symmetric  system,  we  observe 
that  the  optimum  power  assignment  scheme  gains  more  power 
savings  in  the  asymmetric  system  than  in  the  symmetric  system. 


The  worst  SINR  vs.  variable  x,  K=2,  R  .  =30dB 

total 


Fig.  6.  Maximization  of  the  minimum  SINR  among  all  source-destination  pairs 
under  given  total  power  budget  Ft  otai  =  30  dB  (  K  —  2  ). 


TABLE  IV 

Optimum  Power  Assignment  that  Maximizes  the  Minimum  SINR 
Among  Two  Source-Destination  Pairs  Under  Given  Total  Power 
Budget  FtotaI  =  30  dB,  by  Using  the  Proposed  Method  and  the 
Exhaustive  Search  Method  [K  =  2  ) 


Pi(dB) 

P2(dB) 

Pr(dB) 

Asymmetric 

System 

Proposed  method 

19.22 

25.34 

27.59 

Exhaustive  search 

19.37 

25.45 

27.66 

Symmetric 

System 

Proposed  method 

24.49 

24.49 

26.41 

Exhaustive  search 

24.64 

24.64 

26.50 

Also,  we  observe  that  the  total  power  consumption  of  the  equal 
power  assignment  scheme  is  the  same  in  each  of  the  asymmetric 
and  symmetric  systems,  which  is  consistent  with  our  previous 
discussion  that  the  equal  power  assignment  depends  only  on  the 
most  challenging/weakest  source-destination  pair. 

B.  Max-Min  SINR  Optimization 

In  the  second  set  of  numerical  studies,  we  illustrate  the  op¬ 
timum  power  assignment  algorithm  that  maximizes  the  min¬ 
imum  SINR  among  all  source-destination  pairs  under  a  given 
total  power  budget.  We  consider  a  system  with  two  source-des¬ 
tination  pairs  and  one  relay,  i.e.,  K  —  2,  with  a  total  power 
budget  Ptotal  =  30  dB.  We  assume  that  the  cross-correlation  of 
the  two  source  codes  is  p  =  0.25.  We  consider  both  an  asym¬ 
metric  case  (the  system  setup  is  the  same  as  that  in  Fig.  2)  and  a 
symmetric  case  (the  system  setup  is  the  same  as  that  in  Fig.  3). 

Fig.  6  shows  the  maximization  of  the  minimum  SINR  with 
varying  parameter  x  >  0  for  both  the  asymmetric  and  sym¬ 
metric  systems,  respectively.  The  optimal  values  of  the  param¬ 
eter  that  maximize  the  minimum  SINR  of  all  source-destination 
pairs  are  x*  =  1.16  for  the  asymmetric  system  and  x*  =  0.69 
for  the  symmetric  system.  For  both  the  asymmetric  and  sym¬ 
metric  cases,  the  optimum  power  assignments  Pi,  P2  and  Pr 
are  determined  based  on  (58)  and  (59)  with  the  corresponding 
optimal  value  x*,  listed  in  Table  IV.  In  the  table,  we  also  com¬ 
pare  the  optimal  power  values  obtained  by  the  proposed  approx¬ 
imation  method  and  those  obtained  by  exhaustive  search  based 
on  the  optimization  in  (51).  We  observe  that  the  optimal  power 
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The  max-min  SINR  vs.  total  power  budget,  K=2,  non-orthogonal 


Fig.  7.  Comparison  of  the  minimum  SINR  resulting  from  the  proposed  op¬ 
timum  power  assignment  scheme  and  the  equal  power  assignment  scheme  with 
any  given  total  power  budget  ( Ii  —  2  ). 


values  of  the  power  assignment  obtained  by  the  two  methods 
are  indistinguishable.  We  note  that  for  the  symmetric  case,  the 
sources  are  allocated  equal  power  under  the  max-min  SINR  op¬ 
timization  scheme,  while  the  relay  utilizes  transmission  power 
which  is  some  what  different  from  that  allocated  to  the  sources. 

Finally,  in  Fig.  7,  we  compare  the  minimum  SINR  resulting 
from  the  proposed  optimum  power  assignment  scheme  and  the 
minimum  SINR  resulting  from  the  equal  power  assignment 
scheme.  We  consider  both  an  asymmetric  case  and  a  symmetric 
case,  and  the  system  setup  is  the  same  as  that  in  Fig.  6.  We  can 
see  that  in  the  asymmetric  case,  the  minimum  SINR  of  the  two 
source-destination  pairs  is  improved  significantly  (by  at  least 
2  dB)  when  we  use  the  proposed  optimum  power  assignment 
scheme  instead  of  the  equal  power  assignment  scheme.  In 
the  symmetric  case,  we  can  see  that  the  performance  of  the 
equal  power  assignment  scheme  is  very  close  to  that  of  the 
optimum  power  assignment  scheme  (actually,  in  this  case  all 
the  source-destination  pairs  have  the  same  performance). 

VII.  Conclusion 

In  this  paper,  we  analyzed  and  optimized  a  multi- source 
multi-destination  relay  network  where  a  relay  amplifies  and 
forwards  simultaneously  the  signals  received  from  all  sources. 
We  developed  optimum  power  assignment  schemes  for  two 
scenarios.  The  first  scheme  minimizes  the  total  power  consump¬ 
tion  of  all  sources  and  the  relay  under  the  constraint  that  the 
SINR  requirement  of  each  source-destination  pair  is  satisfied, 
while  the  second  scheme  maximizes  the  minimum  SINR  of 
all  source-destination  pairs  with  any  given  total  power  budget. 
Clearly,  both  optimization  problems  as  stated  above  involve  K 
power  variables,  where  K  is  the  number  of  source-destination 
pairs  in  the  network,  which  implies  that  an  exhaustive  search 
approach  is  prohibitive  for  large  K.  In  this  paper,  we  derived 
an  asymptotically  tight  approximation  of  the  SINR  that  allows 
us  to  reformulate  the  original  optimization  problems,  and  even¬ 
tually  reduce  them  to  single-parameter  optimization  problems, 


which  can  be  easily  solved  by  numerical  search  of  the  single  pa¬ 
rameter.  Then,  the  corresponding  optimal  transmission  power 
at  each  source  and  at  the  relay  can  be  calculated  directly.  Each 
of  the  proposed  optimization  scheme  is  scalable  and  the  power 
assignment  algorithm  has  the  same  optimization  complexity 
for  any  number  of  source-destination  pairs  in  the  network. 
For  the  special  case  of  transmission  over  orthogonal  channels 
under  quasi- static  channel  conditions,  we  were  able  to  further 
simplify  the  proposed  single-variable  optimization  scheme,  and 
obtain  an  analytical  solution  for  a  symmetric  system.  Numerical 
studies  were  provided  to  illustrate  and  validate  our  theoretical 
developments. 
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